Aperiodic conductivity oscillations in quasi-ballistic graphene heterojunctions 
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We observe conductivity oscillations with aperiodic spacing to only one side of the tunneling 
current in a dual-gated graphene field effect transistor with an n-p-n type potential barrier. The 
spacing and width of these oscillatoins were found to be inconsistent with pure Farbry-Perot-type 
interferences, but are in quantitative agreement with theoretical predictions that attribute them to 
resonant tunneling through quasi-bound impurity states. This observation may be understood as 
another signature of Klein tunneling in graphene heterojunctions and is of importance for future 
development and modeling of graphene based nanoelectronic devices. 



Graphene is a two-dimensional monolayer of carbon 
atoms that results in a zero-gap semiconductor with out- 
standing electronic [l| and thermal properties Q. Un- 
like in most conventional semiconductors, charge carri- 
ers in graphene obey the Dirac equation and are capable 
of ballistic and coherent transport 0, [|[ , and Vese- 
lago lensing [9(. One of the most promising devices for 
applications in nanoelectronics is the graphene field ef- 
fect transistor (GFET), which was shown to be capable 
of ultra high frequency (100 GHz) operation (lfj| . 

Locally gated GFETs give rise to more complex ar- 
chitectures such as n-p-n heterojunctions, which can be 
realized without physically doping the underlying mate- 
rial. In the conventional transport regime the tunneling 
current across a potential barrier decreases with increas- 
ing barrier energy. In contrast, the chiral Fermions in 
graphene have been predicted to tunnel through poten- 
tial barriers with near unitary probability [1, [ll|. As a 
result, an increasing tunneling current with increasing 
barrier energy is expected, analogous to Klein tunneling 
in quantum electrodynamics. Furthermore, top gated 
graphene heterojunctions form a Fabry-Perot (FP) type 
cavity for electron waves, which undergo oscillations in 
the top gated region due to multiple reflections within 
the barrier. A 7r-phase shift of the FP oscillations in a 
magnetic field was recently observed and understood to 
be the signature of Klein tunneling in graphene hetero- 
junctions 0,0, EH- Recently, Rossi et al. predicted that 
the residual impurity concentration in partially disor- 
dered junctions gives rise to a non-negligible scattering 
potential V sc , which causes broad and aperiodic con- 
ductivity oscillations, which might be understood as an- 
other signature of Klein tunneling (l3j . 

Here we report the experimental observation of aperi- 
odic conductivity oscillations in a quasi-ballistic GFET. 
We analyze the spacing and the width of these oscilla- 
tions, and find that both are better explained by the res- 
onant tunneling model, and cannot be attributed solely 
to Fabry-Perot oscillations. We further analyze the dis- 
order potential introduced by charged impurities and 
phonons at elevated temperatures. 
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Figure 1: (a) Schematic of device geometry and contacting 
scheme, (b) Potential profile created by the global backgate 
and the local top gate as a function of the device length 
(transport direction) creating an n-p-n junction, (c) Room 
temperature (black line) and cryogenic (blue line) conductiv- 
ities of the device. Inset: Optical micrograph of the GFET 
- the graphene flake has been outlined for clarity. 



A schematic of the dual-gated GFET and the electri- 
cal biasing scheme are shown in Fig la. 

The graphene flakes were prepared by micromechan- 
ical exfoliation of natural graphite onto a p++ Si wafer 
with a thermally grown 300 nm Si02 dielectric. The 
flakes were identified as being sin gle layer from their 
characteristic Raman spectrum (lil Il5j . Electron beam 
lithography was utilized to pattern the graphene flake 
and to define electrical Cr/Au contacts. A 10 nm AI2O3 
layer was then evaporated to provide a gate oxide for the 
subsequent deposition of a 100 nm wide top-gate. An 
optical micrograph of a device similar to the one used 
in these experiments is shown in the inset to Fig. lc. 
The top gate was used to apply a local electrostatic po- 
tential, thereby creating an n-p-n junction as shown in 
Fig. lb, with a partially graded junction of width D u . 

The visibility of conductivity oscillations depends 
strongly on the underlying mobility, coherence length, 
and the disorder potential. We first characterize the 
relevant transport parameters in our GFET device by 
grounding the top gate to the drain electrode. Conse- 
quently, carrier transport in the resultant structure is 
similar to transport in a graphene nanoribbon (GNR). 
A constant 100 mV source drain bias V s d was applied 
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while the global backgate bias Vb g was varied from -40 V 
to +40 V. The relatively high backgate bias was re- 
quired since a 300 nm thick dielectric is necessary to 
provide good optical contrast for the purpose of locat- 
ing graphene flakes. The mobility can be estimated 
from the data in Fig.lc using fj, = (enp) , where 
n = C ox {Vbg - Vmrac)/e, V Dirac the voltage at the 
charge neutrality point, and C ox = 115 aF///m 2 the 
oxide capacitance [H, [ljl [l7j|. At carrier densities of 
1 — 2 x 10 12 cm~ 2 , corresponding to back gate voltages 
of 30-40 V, we estimate a room temperature (RT) mobil- 
ity of 1120 cm 2 V _1 s _1 and a cryogenic (5 K) mobility 
of 3300 cm 2 V _1 s~ 1 . The corresponding ballistic mean 
free paths of l e = 50 nm at RT and l e = 110 nm at 5K 
were estimated from the scattering time [l6j]. Note that 
the a — Vbg curve in Fig. lc shows some kinks that most 
likely originate from tunneling through trapped states 
which originate from Fermi-level pinning of the local 
potential at the impurities. These kinks occur on both 
sides of the global conductivity minimum, and their po- 
sition changes with each cooldown event. 

Furthermore, the 2D sheet carrier density An s was es- 
timated from the geometry of the device and known ma- 
terial properties using An s = e e r (V bg - V Dlrac )/ed ox , 
where d ox = 300 nm is the oxide thickness, e r is the di- 
electric constant of Si02- Similarly, the corresponding 
top-gate carrier density Autg can be calculated from 
the sheet carrier density An s [![. 

Unlike top and back gate potentials, An s provides 
unambiguous information about the underlying trans- 
port since An s is zero at the charge neutrality point, 
which is not necessarily located at Vb g = 0. The lo- 
cation of the charge neutrality point away from Vb g = 
originates from the presence of charged impurities, 
which contribute to the conductivity a C i according to 
[H, E3: cr C i (n) = C lmp \An s /n lmp \ , where C imp = 
5 x 10 15 V~ 1 s~ 1 is a constant related to the screened 
Coulomb potential [20j . and nj mp is the impurity den- 
sity. From this equation we find rii mp — 6 x 10 11 cm -2 
in our device, the knowledge of which is crucial in iden- 
tifying the transport regime, and correlating it to the 
visibility of conductivity oscillations. We further calcu- 
late the ft parameter, given by /3 = n / n i 3 ^ 2 , where n' is 
the slope of the density profile around zero density, and 
rii = e/[ih, which differentiates between diffusive and 
ballistic regimes [2l|. Values of /3 -C 1 indicate purely 
diffusive transport, whereas values of /3 3> 1 are char- 
acteristic of the ballistic regime. In our device (3 = 3.7 
at 5 K, which is indicative of the quasi-ballistic regime 
were both ballistic and diffusive transport contributes 
to the conductivity. 

Following the initial characterization of the GNR, the 
main results have been achieved by applying a bias Vtg 
to the top gate thereby creating the electrostatic poten- 
tial shown in Fig. lb. Using this configuration, we ob- 
served an increasing tunneling current with increasing 
barrier height in the vicinity of the charge neutrality 
point, as shown in Fig. 2a. In addition, the onset of 
up to three conductivity oscillation minima (resistivity 
maxima) are visible in an aperiodic spacing to only one 
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Figure 2: (a) Conductivity as a function of back and top 
gate bias showing conductivity oscillations characteristic of 
Klein tunneling. The lines mark conductivity minima, (b) 
Positions of the peak minima as a function of top and back 
gate bias. Data recorded at 5 K. 



side of the global conductivity minimum An tg = 0. 

Similar conductance oscillations have been previously 
observed 0, [It], HH, and analyzed in the context of 
FP oscillations. The oscillations observed in our de- 
vice cannot be explained solely by the FP model. In 
the FP model, where the k-vector is affected by the 
geometric boundary, the magnitude of the spacing be- 
tween successive oscillations S%j can be approximated 
by the condition kp (2L) = 2tt. Consequently the 



peak spacing becomes Sij — 



w 

2d, and is thus constant 
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on the order of 5 x 10 cm 



/nno/Lc, where L c = 
[22], which would be 
in our device. How- 



ever, in the presence of a scattering potential V sc due 
to impurity states, the phase shift of the interference 
&wkb is given by 9 W kb = - J V sc (x' ,y) dx' , where 
the scattering potential V sc is proportional to V sc ~ 

ill 



V d {v) + V tg (r) + £ /dV 



|r-r'| 



[13J . Using the above 



expression for the impurity potential Rossi et al. ob- 
tain peak spacings of 5\2 = 0.85xl0 12 cm " and 823 = 
1.0 xlO 12 cm at an impurity concentration ni mp = 
5 x 10 11 cm- 2 [13]. In our device S 12 = 0.8 x 10 12 cm~ 2 
and 623 — 1.1 x 10 
concentration of n,;,. 
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cm at the estimated impurity 
i = 6x 10 11 cnr 2 . 

Furthermore, the observed aperiodic spacing, i.e., 
812 7^ <^23 as shown in Fig. 2b, can be accounted for 
by the fact that V sc is a function of topgate bias as well 
as the carrier concentration inside the junction. This 
feature cannot be explained solely by the FP model but 
it is present in the self-consistent simulations using the 
above expression for the scattering potential. In addi- 
tion, the width of our oscillations is larger than the the- 
oretically predicted width of FP oscillations [l2j, and 
the experimentally measured width in cleaner devices 
[E, 0] • In contrast, in the presence of V sc , the visibility 
of these oscillations is strongest at low impurity concen- 
trations and decreases with increasing impurity density, 
but are predicted to be still visible at our experimental 
values of n imp = 6x 10 11 cm~ 2 , while both the theoret- 
ical and experimental width are estimated to be about 
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Figure 3: (a) Temperature dependence of the conductivity 
oscillations, (b) Amplitude AI of the first oscillation with 
respect to the minimum between the first and second peak, 
(c) Temperature dependence of the resistivity as a function 
of several back-gate biases. 

0.4 x 10 12 cm -2 , and thus comparable. Consequently, 
the magnitude of oscillation spacing and width in our 
device is best explained by resonant tunneling through 
quasi-bound impurities. 

Finally, we analyze the temperature dependence of 
the conductivity oscillations to study the influence of 
the degradation of the ballistic mean free path and polar 
optical phonons on the observed oscillations. The mag- 
nitude of oscillations diminishes at higher temperatures 
as shown in Fig. 3a. Figure 3b shows a log-log plot of 
the amplitude of the first resistance oscillation peak AI 
as a function of temperature, which vanishes at about 
85 K. To understand what happens at this particular 



temperature we recorded the temperature dependence 
of the GFET resistivity with the top-gate grounded as 
shown in Fig. 3c. The observed exponential increase 
in resistivity is consistent with previous investigations 
fill . [23l | . The fitted solid curves in Fig. 3c corre- 
spond to the phonon contribution to the resistivity p prl 
as given by p ph = p (Vbg) + Pla (T) +ppo (Vbg,T), 
where po (Vbg) is the residual resistivity, and p^A and 
Ppo are due to acoustic and polar optical phonons, re- 
spectively (lq j. The only free parameter in p^A is the 
acoustic deformation potential Da, which we extract 
from the linear part in Fig. 3c to be Da = 15 ± 3 eV. 
The contribution to the total phonon resistivity due to 
Ppo was analyzed following Ref. (24|, and is in good 
agreement with our data (see fit in Fig 3c). 

Thus, the underlying scattering mechanism most 
likely originates from polar phonon injection from the 
underlying SiC>2, which set in above 85 K. These 
phonons cause carrier scattering and thereby degrade 
the ballistic mean free path. As a consequence, trans- 
port through the barrier becomes purely diffusive, since 
l e = 110 nm at 5 K degraded to values smaller than 
the top gate length of 100 nm which defines the n-p-n 
junction. Therefore, the conductance oscillations vanish 
at about 85 K when the GFET transitions from quasi 
ballistic to diffusive transport through the barrier due 
to the onset of polar optical phonons. 

In summary, we fabricated an GFET with an n-p-n 
type potential barrier and observed aperiodic conduc- 
tivity oscillations in the quasi-ballistic regime (/3 = 3.7). 
The peak spacing cannot be explained solely by the FP 
model, but is correctly predicted when resonant tun- 
neling through impurity states is taken into account, in 
agreement with recent theoretical predictions. The ob- 
servation of resonant tunneling through impurity states 
may be understood as another signature of Klein tun- 
neling in graphene hetero junctions and is of importance 
for future development of high performance GFETs. 
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